A Software Package for the Numerical Integration of ODEs by Means of High-Order Taylor Methods
نویسندگان
چکیده
This paper revisits the Taylor method for the numerical integration of initial value problems of Ordinary Differential Equations (ODEs). The main goal is to show that the Taylor method can be competitive, both in speed and accuracy, with the standard methods. To this end, we present a computer program that outputs an specific numerical integrator for a given set of ODEs. The generated code includes adaptive selection of order and step size at run time. The package provides support for several extended precision arithmetics, including user-defined types. The paper discusses the performance of the resulting integrator in some examples, showing that it is a very competitive method in many situations. This is specially true for integrations that require extended precision arithmetic. The main drawback is that the Taylor method is an explicit method, so it has all the limitations of these kind of schemes. For instance, it is not suitable for stiff systems.
منابع مشابه
A software package for the numerical integration of ODE by means of high-order Taylor methods
متن کامل
Numerical solution and simulation of random differential equations with Wiener and compound Poisson Processes
Ordinary differential equations(ODEs) with stochastic processes in their vector field, have lots of applications in science and engineering. The main purpose of this article is to investigate the numerical methods for ODEs with Wiener and Compound Poisson processes in more than one dimension. Ordinary differential equations with Ito diffusion which is a solution of an Ito stochastic differentia...
متن کاملChebyshev Spectral Collocation Method for Computing Numerical Solution of Telegraph Equation
In this paper, the Chebyshev spectral collocation method(CSCM) for one-dimensional linear hyperbolic telegraph equation is presented. Chebyshev spectral collocation method have become very useful in providing highly accurate solutions to partial differential equations. A straightforward implementation of these methods involves the use of spectral differentiation matrices. Firstly, we transform ...
متن کاملNordsieck representation of high order predictor-corrector Obreshkov methods and their implementation
Predictor-corrector (PC) methods for the numerical solution of stiff ODEs can be extended to include the second derivative of the solution. In this paper, we consider second derivative PC methods with the three-step second derivative Adams-Bashforth as predictor and two-step second derivative Adams-Moulton as corrector which both methods have order six. Implementation of the proposed PC method ...
متن کاملA Novel Finite Difference Method of Order Three for the Third Order Boundary Value Problem in ODEs
In this article we have developed third order exact finite difference method for the numerical solution of third order boundary value problems. We constructed our numerical technique without change in structure of the coefficient matrix of the second-order method in cite{Pand}. We have discussed convergence of the proposed method. Numerical experiments on model test problems approves the simply...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Experimental Mathematics
دوره 14 شماره
صفحات -
تاریخ انتشار 2005